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Cn ' Abstract: Massless interacting scalar fields in de Sitter space have long been known to expe- 

rience large fluctuations over length scales larger than Hubble distances. A similar situation 
'^ ■ arises in condensed matter physics in the vicinity of a critical point, and in this better- 

^r^ • understood situation these large fluctuations indicate the failure in this regime of mean-field 

ly-^ . methods. We argue that for non-Goldstone scalars in de Sitter space, these fluctuations can 

^^ i also be interpreted as signaling the complete breakdown of the semi-classical methods widely 

used throughout cosmology. By power-counting the infrared properties of Feynman graphs 
in de Sitter space we find that for a massive scalar interacting through a A(^^ interaction, 
L^ ■ control over the loop approximation is lost for masses smaller than m ~ y/XH/27:, where H 

^ • is the Hubble scale. We briefly discuss some potential implications for inflationary cosmology. 
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1. Introduction 

Some aspects of quantum fields on de Sitter space remain controversial, long after their first 
investigation more than 30 years ago ||l|, but their potential relevance as an explanation of 
the detailed properties of the fluctuations observed by precision measurements [||, |3| in the 
Cosmic Microwave Background (CMB) radiation has stimulated a recent re-examination of 
the issues |^, |5|, |6|, |^, |^, ^, 10, 11, |l2|, Q (see also [15| and references therein) 



The main interpretational difficulties for de Sitter space arise for massless fields, or for 
those that are very light compared with the Hubble scale: m <^ H. There are two related 
ingredients that complicate calculations with such fields: the presence of various types of 
infrared singularities, and the presence of large fluctuations over extra-Hubble distances. 

For example, if a massless scalar field is prepared in an initial state for which fluctua- 
tions are small, then quantities like ((/)^(t)) grow linearly with cosmic time, t. For massive 
scalars this growth eventually saturates at a t-independent value {(jP') oc H^/m?, which is 



parametrically large if m <C -ff [16|. Since these fluctuations are uncorrelated, {4>{x)4>{y)) ~ 
{(p{x)) {(p{y)) , on scales longer than H~^ and since gradients quickly redshift away, one finds a 
picture in which the field takes an approximately constant value within each Hubble volume, 
with different volumes evolving independently of one another in an essentially random and 



uncorrelated way [17|. 



Although this has long been recognized as the appropriate physical picture, essentially 
all of what we know about fields in de Sitter backgrounds is based on calculations performed 



within the semi-classical approximation. This approximation assumes that classical field the- 
ory captures the dominant physics, with small calculable corrections arising from quantum 
fluctuations. The size of these corrections is believed to be kept small because of their system- 
atic dependence on small dimensionless quantities, like powers of coupling constants, A/(47r)^, 
and (because gravity is non-renormalizable) of small energy ratios, E/Mp |l^. In particular, 
such an approximation underpins calculations of inflaton fluctuations during inflation, and 
their implications for the properties of the CMB at recombination. 

In this paper we argue that for massless scalar fields in de Sitter backgrounds subject 
to non-derivative self-interactions — like l^nt — A c/)^ — the presence of large extra-Hubble 
fluctuations undermines the entire semi-classical approximation. This is because the semi- 
classical approximation is at heart a mean-field approximation, within which a quantum 
field is represented as a dominant classical configuration plus a small quantum fluctuation, 
4){x) = ip{x) + 54>{x). But this kind of a description fundamentally breaks down over distances 
larger than H^^ due to the large fluctuations occurring on these scales. 

Notice that the breakdown of semiclassical methods we have in mind does not merely 
mean that the classical approximation is inadequate, with the situation being saved if we 
compute just a few more loops than usual. Instead, for massless scalars the danger is that the 
assumption that higher loops are suppressed by a small quantity that breaks down, meaning 
that all orders in the semiclassical expansion have similar sizes. This makes semiclassical 
calculations inherently unreliable because the truncation of the loop expansion omits contri- 
butions that are as large as those that are kept. 

Quantum fleld theory is a cruel but fair master, so (as usual) the formalism contains 
within itself the news about the breakdown of semiclassical methods. The messenger is in 
this case the infrared divergences that commonly plague de Sitter calculations with massless 
scalars. These indicate a singular dependence on the mass in the more general case of a 
massive, but very light, scalar field. We argue that at least some of these divergences reflect 
the dominance of fluctuations over the contributions of the classical background, pointing to 
a fundamental breakdown of mean-field methods. 

To make our claim precise, we consider a scalar field in de Sitter space, that self-interacts 
through a quartic scalar potential: V = | mgc/)^ -|- ^Xfj)'^. Perturbing in A and using prop- 
agators for fields of mass mo shows that in the small-mass limit the usual loop-suppressing 
factor of [A/(27r)^] for each loop in an L-loop graph is systematically modified by factors of 
{H'^ /ijiq), indicating that higher loops are not suppressed once the scalar mass is sufficiently 
light. 

We explicitly identify contributions to L-loop correlators that are proportional to 



An'^mQ 



which, taken at face value, would indicate that perturbation theory fails once rriQ ~ vA {H/27r)'^. 
However we argue that the breakdown of perturbative methods at this higher mass arises 
because the physical mass scale that cuts off IR effects is really m? = mi + 6m'^, with 



5m? ~ XH^/tjiq, and so it is only the expansion H"^ /rm? = (i:/^^/mg)[l — 5m? /rriQ + • • •] that 
breaks down when ttiq ~ vA(-ff/27r)^. This particular breakdown can be resummed: that 
is, it can be removed by reorganizing the perturbative expansion so that the unperturbed 
lagrangian involves the mass term m?(f? rather than rniQCJ?. This particular reorganization 
does not also remove the potential breakdown of perturbation theory at rrig ~ X{H/2Tr)'^. 

Finite-temperature field theory provides a well-understood precedent for these conclu- 
sions. The small-momentum limit of the Bose-Einstein distribution function, ng{k) oc T/k, 
implies that infrared divergences are stronger at finite temperature than they are at zero 
temperature. As a result, an L-loop Feynman graph for a self-interacting scalar field at finite 
temperature comes with the systematic factor (AT/47r^?7io) , again indicating a breakdown 
of the loop expansion if mo < AT/(27r)^. In the thermal case it is also known that the loss 
of perturbative control is only partial when ttiq = since then the complete thermal mass 
is m^ = 5m? oc AT^ and the IR divergences can be weakened by reorganizing perturbation 
theory so that the unperturbed fields have this mass. In this case the same power-counting 
leads to the loop factor XT/m oc vA, leading to controlled (but non-analytic) results for small 
A. This same resummation fails, however, if ttiq < is adjusted so that tt? = 0, leading to a 
bona fide breakdown of semiclassical methods, such as is known to occur in condensed matter 
systems in the vicinity of a critical point. The breakdown of the loop expansion in this case 



underlies the well-known failure of mean-field methods to describe critical exponents |19]. 
The physical origin of this breakdown is the dominance of large fluctuations near the critical 
point, similar to the fluctuations found in de Sitter space. 

The analogy between thermal field theory and de Sitter space is robust, apart from one 
potentially important difference: in de Sitter space the difference 5m'^ = m? — uiq oc XH^/uiq 
itself depends singularly on ttzq, unlike for thermal field theory where 5m'^ oc AT^. Because 
of this difference the map between the two cases is not simply to think of the Hubble scale 
as a temperature. In particular, although m.^ can be adjusted to vanish in the thermal case 
by appropriately choosing ttiq, it is not clear that this can be done in the de Sitter situation. 
We discuss in the conclusions the open problem of the extent to which these same regimes of 



resummation also apply in the de Sitter case (however, see pO|| ). 

Although gravity resembles a massless scalar in many ways, we emphasize that we do not 
expect this same failure to arise for pure gravity in a de Sitter background. The difference 
arises because the gravitational self-interactions are largely derivative couplings, and so are 
typically not as divergent in the infrared. Since Goldstone scalars similarly couple derivatively, 
they also need not share the same mean- field breakdown. Massless spin-1 fields can couple 
without derivatives, and infrared effects are also known to ruin mean- field methods for a hot 
plasma of charged particles interacting through gauge interactions |^. We leave it open 
whether a similar breakdown occurs on de Sitter space, but any such a failure would require 
the existence of very light charged degrees of freedom to survive the exponential de Sitter 
red-shifting. There is some work for SQED in de Sitter space that shows that the photon can 



indeed get a mass there (see [22| for a review). 



The remainder of the paper is organized as follows. The next section, §2, briefly reviews 



the power-counting of infrared divergences for self-interacting scalar fields at finite temper- 
ature, to show how these track the breakdown of mean-field methods. §3 then provides a 
similar power-counting for self-interacting scalars in de Sitter space. Dimensional analysis is 
first used to argue that the worst divergences are logarithmic for any correlation functions. 
For a massive propagator, this translates into inverse powers of mass at each loop order. A 
class of graphs is then displayed that verifies this dependence through explicit calculation. 
Finally, §4 draws some preliminary implications for inflationary calculations, and summarizes 
our conclusions. 

2. IR divergences at finite temperature 

We first review the breakdown of mean-field methods near a critical point. To this end we 
estimate the size of the contribution to physical amplitudes of the infrared divergences that 
arise due to long- wavelength fluctuations. The purpose is to sketch why these long- wavelength 
fluctuations eventually dominate the suppressions due to small couplings at any order in a 
perturbative expansion. 

Consider then a massive scalar field with a quartic self interaction in flat spacetime 

C = -'-d,<pd^4>-^^'-^<t>\ (2.1) 

heated to a nonzero temperature T. In Euclidean signature the scalar propagator is 

i 1 

G'n(p) = -JTT^ -^—. 2~: 2 ' (^•^) 

(27r)'^ Pn + P + rriQ 

where Pn = 27rn//3 = InnT, with n = 0, ±1, ±2, • • •, ensures periodicity in imaginary time 
r-^ r-h/3. 

Consider now a Feynman graph that involves E external lines, / internal lines and V 
vertices, and (to start) suppose there are no insertions of the mass counter-term vertex 5m?. 
Then the parameters E, I and V are related by the identity (conservation of ends) E+2I = AV 
and the definition of the number of loops, L = I — V + 1. Using these to eliminate / and V 
then gives 

2V = E + 2{L-l) and 21 = E + A{L - I) . (2.3) 

If /c ~ vfiQ denotes a typical external momentum scale, the resulting amplitude is schemat- 
ically given by 



^°(k,T)~ 



1 ^^ 
A(2^)3<53(p + k)-5„„, 
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(2.4) 



All of the delta functions perform one of the integrals, except for one which expresses the 
overall momentum conservation delta function. 



The most IR singular part of the result corresponds to the n = term of the sums, so if 
the IR convergence occurs at scales p ~ tuq (more about this below) then for k ~ niQ <^ T 
the most singular dependence on rriQ has the form 
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(27r)2mo 



(2.5) 



This is the main result, expressing the most singular small-mo limit of an L-loop contribution. 

Perhaps the biggest surprise in the estimate (|2.5| ) is the systematic appearance (when 
mo <C T) of the factor T/mo at each loop order. Once m-Q ^ A^T^/lGvr^ this eventually 
undermines the semiclassical loop expansion itself, whose validity ultimately rests on the 
suppression of loops by powers of A <C 1. 

What happens in the limit uiq — )• can be understood by reorganizing perturbation 
theory to recognize that finite-temperature effects also contribute to the scalar mass, and so 
can themselves suppress the total failure of an expansion in powers of A. This reorganization 
can be made explicit by rewriting the lagrangian, eq. (p.l|), as 
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(2.6) 



where rr? = m^ + 5m? is the physical scale responsible for saturating IR fluctuations, and 
5rn? is a perturbatively small, but calculable, thermal mass shift, whose leading form is 
5rn? ~ cXT"^ /A.'k'^ for c a positive constant of order unity. 

Now consider repeating the above power-counting argument, but also inserting n factors 
of the mass counter-term, 5m? , into the result. Since each such insertion also adds a new 
propagator, it contributes an amount ~ 5m? /'p?, so the same dimensional argument as given 
above for the leading IR behaviour results in the revised estimate 



A^,{k,T)c^5\\.) 
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(27r)2m 
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m^ 



(2.7) 



In the limit where the zero-temperature mass vanishes, mg = 0, then m? = 5m? ~ cXT'^ /Att"^ , 
and so all of the factors 5m? jm? are order one. But the point of this reorganization is that 
these new order-unity interactions proportional to 5m? systematically cancel order-unity parts 
of the self-energy corrections in Feynman graphs, leaving a result that is suppressed by a net 
power of A. 



For instance, the one-particle reducible contribution to the 2-point function at zero ex- 
ternal momentum at two loops has an IR singularity proportional to }?T'^/m?. But this 
precisely cancels the contributions of the one-particle reducible graphs involving Sm? inser- 
tions, to leave a residual contribution that is of order /(^T"^ /m oc A^' ^T^, in agreement with 
standard calculations. In general, using Sm? = m?' ~ A (T/27r)^ in eq. (2.7) shows that the 
small parameter that suppresses each loop in this case is 

AT \^ ( V\\^ , . 



iTxYm) \ 27r 

The contribution of higher- loop graphs remain small for small A, although each additional loop 
now costs a factor of vA rather than A itself. This initially surprising non-analytic dependence 
on the coupling A can arise within perturbation theory because of the reorganization of 
perturbation theory that is implicit in using the loop-generated mass to cut off the infrared 
divergences. For scalars it is known that the shift in mass completely resolves the apparent 
breakdown of perturbation theory pl[| , though the same is not true for massless gauge bosons 
in hot plasmas. 

A situation where the above resummation does not salvage the semiclassical loop ex- 
pansion is when the zero-temperature mass is negative and adjusted in precisely the way 
required to ensure that the full finite-temperature mass vanishes for a particular tempera- 
ture: m-Q = —brv? ~ — cAr^/47r^, so that rv? ~ 0. This is what happens at a critical point, 
T = Tc, such as appears at the termination of a phase-coexistence curve. In such a case the 
vanishing of m in the estimate (|2.7| ) indicates that control over the loop expansion breaks 
down completely, in agreement with the well-known failure of mean-field methods to compute 
critical exponents near critical points (where m — )• 0) ||l^. The root of the problem with the 
loop expansion in this case lies in its mean-field nature, since the large fluctuations allowed by 
the massless fluctuations at a critical point invalidate the expansion about a large background 
field. 

3. IR divergences in dS space 

We next repeat the above power-counting argument for calculations of scalar-field correlations 
in a fixed background de Sitter spacetime. To this end we work within the 'in-in' formalism 
[p^ ], and imagine computing a correlation function for (j){x) for a minimally coupled scalar 
field that self-interacts through the potential 

V{cP) = Vo + ^cp' + ^^^'. (3.1) 

We take the background cosmological constant, Vq, to be sufhciently large as to dominate 
the quantum fluctuations of cp, allowing the background curvature to be regarded as a fixed 
geometry ds^ = a^{T) (— dr^ -|- Sijdx'^dx^^ , where r is conformal time and the de Sitter scale 
factor is a = e^* = —1/{Ht). 



3.1 In-in correlations 

Within the in-in formahsm the path integration is over a duphcate set of field configurations, 
0+ and (f)-, that correspond to the two time-paths to temporal infinity that arise, but it is 
useful to instead use the two combinations 0c = \{4'+ + 0-) and cj)/\ = (j)j^ — cj)^. Since the 
path integral is weighted by the integrand exp [iS{(j)j^) — iS'((/>-)], in terms of the fields (j)c 
and (/>A the scalar self-interaction appearing within the path integral therefore has the form 

= m2</.c(/.A + ^ (4(/.^<Aa + (l)c(t>\) , (3.2) 

so these are the interactions whose vertices we follow in any particular Feynman graph. 

The internal lines of the graph represent the correlations of these fields: {(pc'Pc)-, {4'A<Pc) 
and {4>c4'a) (part of the utility of the combinations 0c and 0a is the vanishing of the au- 
tocorrelation (0a0a) = 0). The correlator {4>c{ti)4'c{t2)) for a massless scalar propagator 
becomes 



G°(/c,Ti,r2) = -TTg |(1 + A;^riT2)cos[/c(Ti -r2)] +/c(ti -r2)sin[/c(ri -T2)]| 
where the last line specializes to the long- wavelength, super-Hubble limit. 



-kT = ^<^l. (3.4) 

an 



The retarded correlator, (0c(''"i)0a('^2))) is similarly 



G%{k,Ti,T2) = 9{ti -r2)-p- |(l + /c^Tir2)sin[A;(Ti - T2)] - A;(ri - T2) cos[A;(ti -r2)]| 

{rf-rl){l + 0[{knf]}. (3.5) 



- 0{ti - T2) 



3 

3.2 Power-counting 

For the purposes of power-counting imagine computing a correlation function involving Nc 
powers of 0c and A'^a powers of 0a: {4>^4>c'^). A generic Feynman graph contributing 
to such a quantity involves Ir internal lines involving retarded (or advanced) propagators, 
{(Pc'Pa) and Ic lines representing {(I)c4>c) propagators; linking Vc vertices describing the 0c0A 
interaction, and Va vertices built from 0c 0a — see Figure (IT]). If we write m-Q = m? — 5w? 
then graphs also involving insertions of 5m? have additional factors, but for the moment we 
ignore these. 

There are also Ec external lines corresponding to the G^ propagators, each of which 
must attach to one of the Nc (t>c fields in the correlation. Due to the off-diagonal nature of 
the external retarded propagators these can connect to either an external 0c or 0a field. If 




)AN 



0c Va -X 



Vc Va 4'^ 4'^ 

Figure 1: Labelling for vertices and external lines. For external lines the X indicates which end 
attaches to the vertex. 

there are Ej^ such external propagators we denote by E}^ (respectively E}^ ) the number of 
such propagators that connect to an external (pc (respectively 0a) field, as shown in Figure 
(1). Clearly ^i^^ +Ec = Nc and e'^'^ = Na, and so e'^^ + e'^^ +Ec = Ei, + Ec = Nc + Na- 
The numbers of propagators and vertices are related by the definition of the number of 
loops 

L-{Ic + In) + {Vc +Va) = 1 (3.6) 

and the two conditions expressing the conservation of 'C and 'A' type ends, 

Ec + 2Ic + Iu + E'i^ = Wc + Va 

In + 4'^ = Vc + SFa . (3.7) 

These last two use the fact that 'i?' type propagators (involving retarded or advanced Greens 
functions) are mixed correlations that have one (pc end and one (pA end, while 'C type 
propagators involve (pc fields at both ends. It is useful to use these three identities to eliminate 
Vc, Ir and Ic, to get 

Ir = 1 [e^ + Ec- 4'^) +L-1 + 2VA 

Ic = E^^^ +L-1- 2yA (3.8) 



1 



p(2) ^ F ^ f(1) 



Vc = - [E'k' + Ec + E);>) +L-1-VA. 

Now in a generic Feynman graph, each vertex contributes a factor that has the schematic 
form 

YeTte^c^xJ^{2nfsHp) 

(3.9) 

where r is the time that labels the vertex and the delta function expresses conservation of 
the (co-moving) momenta, p, that meet at the vertex. Similarly, for massless scalars each 
internal 'C type propagator contributes in the IR limit 

W^ —3- > (3-10) 



while each internal 'i?' type propagator gives 

—^ H'^ein - Tj) [rf - Tj) . (3.11) 

Combining results, evaluating a Feynman graph gives a result proportional to 

A{ki,Ti; . . . ■kE,TE) - {2TTf5^{k)Aiki,Ti; . . . ■,kE,TE) , (3.12) 

where the delta function schematically represents the overall conservation of spatial momen- 
tum. The expression A is given schematically by 



/ fT2\Ec [■ r A^ 1 Vc+Va 



(27r)3 2p3 



[(27r)353(^)]Vc+yA-i(3_-^3) 



d'p e{Tj-Ti) , 3 

(2^ 3 ^^^~^^> 



where the /c's generically stand for external momenta, the p's for loop momenta, and we use 
the early part of alphabet (a, b, c) for external time and the middle part i, j, I for internal times 
associated with vertices. In the Appendix we examine several Feynman graphs explicitly, to 
examine the behaviour of the integrations over the vertex times, Tj, and loop momenta, p. 

3.3 IR behaviour 

We next estimate the most IR-sensitive part of such a de Sitter-space graph. For technical 
simplicity, even though we imagine the IR regulator to come from nonzero m?, we perform 
our estimate by working with rn-Q = and cutting off the divergence at a physical scale Aj^, 
relating this scale to niQ and H in a second step. 

The dependence of the final result. A, on both momenta, k, and external times, Ta, 
complicates the use of dimensional analysis in identifying the dominant infrared divergences 
that can arise in general amplitudes. To make progress we assume the fields being correlated 
are evaluated in position space, in which case the position-space amplitudes are obtained from 
the above expressions by multiplying by an appropriate factor of (27r)^^ /d ^ e . Power- 
counting then simplifies in the special case that all external fields are evaluated at the same 
position, since then de Sitter invariance of the vacuum implies {4>c'^ {x)(j) ^^ {x)) is independent 
of x^. In this case the result cannot depend on k or Ta, and so the dominant divergence must 
be a function only of H and the IR cutoff, A^r. Notice that this argument requires A^ to be a 
time-independent cutoff on physical momenta, since a time-independent cutoff on co-moving 
momentum would introduce a spurious time-dependence into the correlation function. (In 



this respect we follow |lj], and at least for the UV [|14[, but differ from the practice used by 
other recent workers.) 



Given the result, eq. ( |3.13 ), of the last section (and after UV divergences are renormalized 



by including the appropriate counterterm graphs) the dominant IR divergence of a position- 



space amplitude of the above type can be determined by dimensional analysis 
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(3.14) 



with 
P={Nc + Ni^)-Ec-l + {Ic + Ir) -Ic- {Vc + ^a - 1) + {Vc + Va) - ^« - /« = . (3.15) 

Since P = 0, we find that in general every diagram contributing to this type of correlation 
function is at worst log divergent in the infrared. Of course, the above dimensional argument 
cannot in itself distinguish a divergence like InAj^ from (InAj^j)^, and so a more precise 
determination (given below) of the nature of the divergence requires a more detailed estimate. 

3.4 Massive Propagators 

To better parse how this divergence is regulated as a function of the scalar mass, we step 
back and use the small-A; behaviour that is relevant to the massive scalar propagator on de 
Sitter space. In the limit — /cr = k/{all) <^ 1 this becomes 



Gc{k,Tl,T2) 
GR{k,Tl,T2) 



2A;3 



ik'^TlT2) 
i/2 



e{n-T2) — {Tt~'°T^°-Tl^W^ 



(3.16) 
(3.17) 



where eo = tjIq/SH'^. Even though these expressions differ only weakly for small k from the 
massless case, they suffice to cure the IR divergences encountered previously because these 
divergences are only logarithmic. 

The main change that this introduces relative to the above estimates is the conversion of 
the internal-line factor 
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Figure 2: A class of graphs contributing terms of order [XH /m^) to Gc{k,T) 

which, using as before a cutoff, Aj^, on physical momenta gives 
" dp fHY^'° ^ [ P 



km P \Pj \^mJ 

This shows that the small-mass limit converts powers of In A^ into powers of H'^/niQ. This 
conversion of a logarithmic divergence into an inverse power of mo (as opposed to powers of 
In mo) arises because for twq <C H^ the difference between the massive and massless expres- 
sions for Gc{k,T) becomes important only once (— A;r)^"^" deviates from unity, which occurs 
for k < ki, with 

aid 

corresponding to a physical IR cutoff of order A^ ~ k^,/a ~ H e~^^ /'^"^o. 

3.5 Powers of logs 

Because power-counting returned the result that graphs diverge like A^ in the infrared with 
P = 0, we know these divergences are logarithmic. What these arguments do not yet show is 
that this logarithmic divergence worsens order-by-order in the loop expansion. To establish 
this we next identify an explicit class of graphs for which additional powers of H/m arise at 
each order in perturbation theory. 

The explicit class of graphs we choose for this purpose are those containing the successive 
chain graphs of Fig. 0. These diagrams are loop corrections to Gc{k, r) = (</>^) and we denote 
by Gq the L-loop chain diagram contribution to {(j)'^). The 1-loop contribution, G^- , is 
evaluated in great detail in |11, 12| so we just quote the end result in momentum space: 

G'^^\k,T) = -A r dTV(T')4°^(A:,T,r')GL°)(fe,T',T)A(r') (3.21) 

J —00 



where 



A(t)= I ^G^^\k,T,r) (3.22) 



11 



is the loop factor. With the appropriate (ultraviolet) mass counterterm and using the massive 
propagator given by Eq. ( |3.16| ), one finds that the loop factor is time independent and equal 



to m 



A(r) 



2eo \2TrJ \h) 



where ^ is some non-IR physical scale. All together, the 1-loop correction to Gq is 



G'a'{k,T,T)) = G)^>{k,T,T: 



6(2vr)2eovl)'"^"(-^") 



(3.23) 



(3.24) 



At L loops the chain diagrams simply factorize, so that the L-loop chain diagram is just 

iL 



.(i) 



G'c'>{k,T,T)xG'^'{k,T,T 



(0), 



X 



(27r)2eo \H 



^y'"ln{-kT) 



(3.25) 



Combining the factors depending on L then shows that each successive loop comes system- 
atically pre-multiplied by a factor of 






mf 



(3.26) 



indicating the breakdown of the loop expansion once uiq is sufficiently small. 



{L)i 



To make contact with our power-counting estimate above we Fourier transform Gq {k, r) 
with respect to k to obtain the corresponding contribution to {(P'^{t)). Keeping in mind that 
Gc {k,T) oc {kT)'^'^° / k^ , the required integral is 



/ 



kMk, 



{kTY'o ln{-kT) 



"'i) 



L+1 



1 + O(eo) 



(3.27) 



and so the ?TT,Q-dependence of the contribution of the L-loop chain graph is^ 



* 47r^m, 



(3.28) 



At face value this indicates the edge of the perturbative domain lies at itiq ~ vA H'^ /2tt, 
where the contribution of fluctuations to the mass begin to compete with uiq. 

Notice, however, that because the dominant contribution to the Fourier transform comes 



from the A; ~ limit of integration, the mass dependence of eq. ( |3.28| ) is exquisitely sensitive 
to the small-/;: form of Gcik,T), which the above graphs shows has an expansion in powers 
of (A/eg) ln(— /cr). Yet we know that for sufficiently small k this expansion breaks down 
because the logarithm systematically competes with the additional power X/cq. A better 
approximation for Gc{k,T) at small k can be found by resumming the leading logarithms to 
obtain a result that is a series in A/eo, without accompanying factors of ln(— fer). 



^We thank Don Marolf and Ian Morrison for helpful correspondence on this point. 
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This resummation can be done ||T^ using the dynamical renornialization group (DRG) 
[jT^ ], by recognizing that because the integrand is a function of fcr, smah k is related to large 
r. The breakdown of the perturbative expansion at small k can therefore alternatively be 
regarded as being a breakdown due to the presence of the secular growth in r; a breakdown 
that the DRG is designed to resum. A better estimate of the small- m, behaviour of {(/)c{t)) 
can therefore be obtained by using the DRG-improved late-time approximation to Gcik,T) 
that also improves its small-Zc asymptotics. 



The result of this resummation |12] is to replace the expansion 



Gc{k,T)=G^^\k,T)\l + 



A 



6(27r)2eo 



iT-^ 



-kT) + 






ln(— A;r' 



with the DRG-improved result 

Gc{k,T)=G)^'ik,T)i-kT) 



.(0), 



2<5 



i + o{ - 



with 



A 



12(27r)22eo 



(3.29) 



(3.30) 



(3.31) 



The key observation is that because Gq {k,T) oc {H^ / k^){—kT)'^'"-\ the small-A; behavior 
of the DRG-resummed contribution to Gc{k,T) is equivalent to what would arise from a 
small shift e^ ^ e = e^ + 6. Equivalently, this corresponds to a shift of the scalar mass 



rrin 



m 



rrif 



+ 5m^ with [||] 



m 



ml + 



3XH^ 

167r2 



nin 



(3.32) 



where e = m? /2>H^ . 

Fourier transforming this more accurate depiction of the small-A; limit in Gc{k,T) to 
obtain (0^(t)) then gives the following integral, 



[r) ^H'l 



__(A:r)2^oc— , 



(3.33) 



instead of eq. (|3.27D . This reproduces the above series in powers of XH^/niQ once expanded 
using 1/e = l/eo — 6/eQ + - • •. This shows that it is the scale m? that cuts off the IR divergences 
in {cI)'^{t)), suggesting the utility of reorganizing the perturbative expansion so that it is the 
mass m? rather than mg that appears in the unperturbed lagrangian. 
The corrections to the DRG-resummed form are of order 



A XH^ 

— oc TT- 



eo 



mS 



(3.34) 



suggesting that the boundary of the reorganized perturbation theory lies at ttt-q ~ XH'^ /Att'^, 
rather than when m'^ ~ yA H'^/Att'^, as was found above. We regard it to be an open 
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question whether the perturbative problems that arise here can themselves be resummed in 
a controllable fashion. The burden on any proponents of resummation is to show that the 
resummed graphs capture all of the leading 1/m behaviour, in the regime of interest. Ex- 
perience with finite-temperature systems argues that although resummation may be possible 
for some regimes (like ttiq ~ 0) this need not imply it can always be done (such as when 
m^ ~ 0). What is not yet clear is whether it is possible in de Sitter space to reach the regime 
m^ ~ 0, since m^ does not pass through zero for any value of uiq, at least within the domain 
of validity for which eq. ( p. 32 ) holds. 



4. Discussion and Conclusions 

The body of this paper argues that an L-loop contribution to a correlation function for a 
scalar field in de Sitter spacetime with A 0^ self-interactions carries a systematic factor of 
[XH"^ / A-K'^rn'o)^ , indicating a fundamental breakdown of semiclassical methods once mo ^ 
vA H/2tt. The origin of the breakdown of perturbative methods is the infrared-singular 
behavior of these graphs which arises due to the large extra-Hubble fluctuations experienced 
by very light scalars in de Sitter space. These fluctuations dominate the contributions to 
correlation functions, invalidating the semiclassical approximation which is at its heart a 
mean-field description. 

This story is qualitatively similar to what happens in finite temperature field theory, 
although the power of coupling constant that defines the boundary of the semiclassical region 
differs. In the finite temperature case, the variance of the field, {(j)^)., goes like T^ while in 
the de Sitter case we find (0^) ~ H^ /m?'. This difference means that the mass, mdym which 
is comparable to the corrections to m scales differently with A in these two cases: for the 
thermal case m\ oc AT^, while for a de Sitter background m\ oc \^''^H'^. It also may 
mean that the structure of de Sitter space precludes access to the regime m?' ~ for any 
choice of ttiq, unlike for finite-temperature systems, see Figure (y). However, as is shown 
above, the two theories have a similar perturbative structure, whose relation is sketched in 
Figure (g). 

A natural question to ask about this perturbative breakdown at small m^ is whether it 
can be resummed (like for ttiq > at finite temperature), or whether it reveals a complete 
breakdown of expansions based on powers of A (like when ttiq < is chosen so that m^ = 
at finite temperature). 

Others have argued that there exist semiclassical methods that capture the leading in- 
frared logs. For secular growing logs, one proposal solves the late-time physics using the 
classical equation of motion llO], in a similar way to the 5N formalism. Another point of 



view uses a stochastic approach to inflation [17], which is argued to capture (and resum 



the leading logs, by generating a dynamical mass [25|. Although such a stochastic approach 
goes beyond mean field, it is not yet clear what combination of small parameters control the 
approximations made in its derivation. 
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Loop size 0{1) 0{\/X) 0{X) 



mc 



'^dyn fT^max 



Thermal AT 



Vat 



de Sitter y/XH A^/^iJ H 



Figure 4: A comparison of the behavior of loop corrections for a scalar field with a quartic interaction 
Xip'^ in a thermal background or de Sitter space as a function of mass. The masses labelled on the 
chart are, from left to right, the mass below which perturbation theory breaks down, the dynamical 
mass generated in the event that the zero temperature or flat space mass was zero, and the maximum 
mass the field can have and still receive thermal/de Sitter fluctuations. 



m 




m 







From the point of view of the ar- 
guments made here, derivative couplings 
are not as dangerous as are those of the 
scalar potential. This is because the mo- 
mentum dependence of these couplings 
tends to ameliorate any IR divergence the 
graph would have otherwise had. This 
means that massless Goldstone bosons on 
de Sitter space would not suffer from the 
same breakdown of perturbation theory 
as does the (j)'^ model considered here. 

We would expect the self-interactions 
of massless gravitons to be similarly be- 



Figure 3: A plot of m^ vs ttiq. The horizontal band 
represents the regime m^ < 0{XH^ /An^) while the 
curve is given by Eq. ( 3.321) . It would be interesting nign so long as these are derivative cou- 
to investigate the case of negative bare mass, which we plings, leading us to expect no pertur- 
do not display here, further bative breakdown for pure gravity on de 

Sitter space. This expectation seems to 
be borne out by ref. [^ and ref. |^ (the latter appeared just as our paper neared comple- 
tion) . These authors study one- loop infrared divergences for pure gravity and gravity coupled 
to scalars in de Sitter and slow-roll spacetimes, but find that IR divergences cancel in the 
absence of the self-interactions of the scalar potential (which appear in their calculations 
as slow-roll parameters). Based on the power-counting arguments presented here we expect 
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this to continue to be true at higher loops, with the scalar self-interactions being the most 
dangerous in the infrared. 

It is interesting to consider in this light what implications our result might have for 
inflationary cosmologies. A complication in directly extracting these for simple single- field 
models is our neglect of classical evolution of the background scalar field and metric, due 
to our use of a simple de Sitter background. Because a homogeneous evolving scalar field 
can be used to define a notion of cosmic time, many of its effects can be gauged away. We 
therefore expect a naive application of the above arguments to simple models are likely to 
cancel from gauge invariant quantities like the curvature fluctuations of physical interest 
for cosmology. Nevertheless, it may be possible to have important infrared effects appear 
in curvature correlations in multi-field models, particularly those involving nontrivial post- 
inflationary dynamics (such as curvaton models). We leave for future work the detailed 
question of whether and how the infrared effects we find above 'propagate' into late-time 
curvature perturbations. 

In the remainder of this section we put these issues aside, and ask what the condition 
> XH^/Att'^ implies for the parameters of a single scalar field described by a quartic 



m 



2 



potential 

V = Vo + ^m''4>^ + ^X(t>''. (4.1) 

For this model we ask how the condition M^ » A(ff/27r)^, where M^ = V" = 710? + ^XcfP', 
compares to other conditions to which inflationary models must be subject. To see what this 
implies for potential inflationary applications, consider two extreme cases: {i) 4> so large that 
y ~ ^ A(/)^; and {ii) 4> small enough that F ~ Vq. (This last example can be regarded as 
part of a model of hybrid inflation |^4| , with inflation ending as another field starts to roll as 
its (/>-dependent squared-mass goes negative.) 

Large- field inflation 

For the large-(/) regime we have F ~ ^ A 0'', and so i:^^ = F/SM^ ~ A 0''/72M2. The slow-roll 
parameters are 

so the edge of the inflationary regime is (psn/Mp ~ 0(1). For (j) larger than this classical 
evolution satisfies 

. V lX<ty^ \/2 



^H VA (AV2V2 Mp 3 



\/A Mp(l) . (4.3) 



Once (j) becomes smaller than H^ fluctuations dominate classical evolution and inflation be- 
comes eternal, which in this case occurs when -^vA Mpcf) < Xcl)^/72M^, or (p^ > (/)|^ ~ 
24\/2 M^/VX. 
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How do the boundaries of the semiclassical approximation compare to these? There 
are two criteria to be satisfied. First, control over the low-energy approximation that un- 
derlies the gravitational loop expansion requires V <C Mp, or cj)'^ <C (1)%^ — '^'^Mp/X. We 
have seen in previous sections that the A loop expansion fails unless M^ ^ XH'^/Att'^ 
or iA(/)2 > X^(P'^/2887r'^M^, and so 0^ < </,|^ ~ {UnMpf/X. Since (/-sfl/Mp ~ 0(1), 
(l)Ei/Mp ~ 0(A-V6), ^^^/Mp ~ 0{X-'^/^) and (pBo/Mp ~ ©(A^VS) ^e have 

(/•Sfl < 0B/ < '?^Hi5 < 4>BD , (4.4) 

and so the condition M^ > X{H/2tt)'^ is parametrically weaker than the condition Mp > V. 
The fact that eternal inflation can occur before a total breakdown of perturbation theory is 
as expected. 

Hybrid inflation 



Consider next a hybrid model |24], involving two scalar fields, (/> and X) interacting through 
the potential 

C/(0, X) = \C{X'- v'f + ^ XV^ + ^ m^^ + 1 ^ '^^ • (4-5) 

In this model the fields start in the trough defined by x = 0, with (p large and rolling towards 
smaller values subject to the effective potential 

y = yQ + i^2^2^^^^4^ (4.6) 

with Vq = 4T C^^- This roll continues until 

<P = ct>sa = ^, (4.7) 

9 

at which point the x mass, /i^ = —C,v^ + g^4P', becomes negative, allowing x fo evolve quickly 
towards the absolute minimum at x = ''^ and = 0. 

The dynamics of the inflaton in this model is governed by the same potential considered 
earlier, but our interest now is in the small-field regime, for which rn^ ^ ■^XcfP' and F ~ Vq 
(and so H^ ~ Vo/3Mp). We assume parameters are chosen to keep (I)sr small enough to 
ensure x remains zero well into this regime. In this case the slow-roll parameters are 

e^-,[-^) and V^^, (4-8) 

so 2e ~ "rf {(j) / Mp)"^ . Clearly inflation only requires m^Mp/Vo <C 1, since the condition 
^ X(fP' <C m? automatically ensures cj) <C Mp/rj. Inflation ends once (j) reaches cpsR = y/C, v/g- 
In this case the validity of the A loop expansion requires m^ 3> A if ^/47r^ ~ AVo/127r^Mp , 

or 

rn^M?; X 
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Notice that this lower bound to the slow-roll parameters is not simply a naturality condition, 
corresponding to a regime for which small loop corrections dominate the smaller classical 
potential. Rather, in this regime the problem is not fixed simply by including one- or two- 
loop corrections to V. Instead it is the entire loop expansion itself that fails. It would be 
of great interest to see whether the semiclassical criterion plays any role in more general 
inflationary contexts. 
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A. Power counting in k-space 

In this appendix we examine the leading infrared behaviour of several some Feynman graphs, 
to follow how the integrations over loop momenta, p, and vertex times, tj, reproduce the 
logarithmic dependence on A^j. We illustrate many of the subtleties that arise when trying 
to do power counting in k space in this context. 

For most diagrams the integral over time either converges well in the infrared {i.e. gives 
terms of order 0(1) -|- C'(— fcr)) or diverges logarithmically, like ln(— /cr) |jl^. Weinberg has 
proven that the log divergence is the worse divergence you can get under a certain set of 
assumptions that are satisfied by our A0^ model [Q]. Nevertheless, for some diagrams the 
time integral give important powers of external momentum which one must keep track of 
to get the power counting right. To see this let us first look at the tree level calculation of 
((/)^). There are two diagrams that contribute to this correlation function at tree level — see 
Fig. (M). For simplicity and illustrative purposes, we take all external momenta and time to 




\ / 




(a) 

Figure 5: The diagram on the left (a) is Ec = 3 Vc = 1 and the diagram on the right (b) has Ec = 1 
and Va — 1- 
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be equal to k and r. The first diagram (diagram (a) with Ec = 3) gives 

D|^(ln(-tT) + 0((-fcTr)) (A.l) 

with n > 0. The integration limits need to be explained; the upper limit arises from the 
9{t — t') in the retarded propagator. By causality, the interaction occurred at any time in the 
past of T. But if the interaction occurs so far in the past that any of the external momenta 
come inside the horizon, —kr' ~ 1, then the corresponding propagator starts oscillating very 
quickly and the integral will in general vanish: f_ ' (• • ■)e~^^'^ ~ 0. The second diagram 
(diagram b) in Fig. (||) on the other hand goes like 

^^Y,mi)+0{(-krr)) (A.2) 

We see that the time integral give the 1/k^ that was naively missing from the second diagram. 
Note that it could be that the 0{1) factors cancel and the diagram would then be ~ ^ such 
that A + B = %. Since we are working with —kr <^ 1, the diagram would be much smaller 
than expected. Our power counting is only good to keep track of the leading piece and it 
assumes no special cancellations. This assumption is completely justified given the absence 
of special symmetry. 

In loop diagrams, we have identified at least three additional complications that arise 
from the integral over internal momentum p and from the entangled time integrals. The first 
subtlety occurs because the momentum of the virtual particle is integrated over all values 
and when p > k, the lower cutoff of the time integral becomes —1/p instead of —1/k. This is 
because the virtual particle stops oscillating at a later time, or leaves the horizon at a later 
time. Coming in from r' ~ — oo, we find that once we reach the time r' ~ — l/fe, the virtual 
propagator is still oscillating and the time integral keeps integrating to until we reach the 
time t' ~ — 1/p- The second subtlety is that the retarded propagator imposes some kind 
of ordering on the internal vertex times tj > r^ > • • • such that the upper bound of some 
integrals are variable of other integrals. Finally the third subtlety arise from the fact that the 
momentum integral may depend on time as well. This could occur for example if one deals 
with infrared effects with a beginning to inflation or if one uses an explicit IR cutoff. 

All these subtleties make the power counting very complicated in k space and this is 
why we limited ourselves to the real space expression in the main text Eq. ( |3.14 ). Still, 



progress can be made to keep track of the leading infrared piece if we make the following set 
of assumptions. We will take all time integral to stop at some external time r and take all 
internal momenta to be smaller than external momenta p < k. We can integrate all internal 
time Tj to external times r^. Finally, the IR physics is a small mass and no time dependence 
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is introduced in the integral over momenta from infrared effects. Applying these simple rules 
to do the power counting of the time integrals in Eq. ( 3.13| ), we get 



time 



dTi 



Vc+Va 



[9{t-t'){ 



/A^^3_^/3)j4''+4'' 



[Ti -T 



iIe 



■irt-r') 



[0{ra-n){Tl 



(2) 



3M^fi 



dTi 



1/k 



Vc+Va 



in) 



3e'-j^^+3Ir~4Vc-4Va 



+ 



(A.3) 



where we have f actor ized the Ej^ part and we have kept only one term (where all r' are Tj) 
typical among many. This gives 



p(2) 1 

time ~ [B{ra - n^r^ - r^)]^- ^^ (0(1) + 0((-fcr)") 



(A.4) 



(2) 

where we combined terms to generate the full correct E^ propagator (it has to) , n > and 
0(1) could be a ln(— fcr) as well. In summary, the main effect of the time integrals is to give 
rise to 1/fc^ '^ and note that the power counting does not keep track of the ln(— fcr). 

To illustrate further the complication of doing these loop diagrams, lets look at the 
following sunset diagram (see Fig P). This is a 2-loop correction to {4>c4'c) which goes like 




Figure 6: The time are from left to right r, ti and T2. The incoming momentum is k. The internal 
momentum from top to bottom are g, p and r = \k — q — p\. 



where we use the notation that ^,-,1 = 0{t — ti). The theta functions give the following 
hierarchy r > ri > T2 and all the propagators are only valid given a series of condition 



-{p,q,k]{T,Ti,T2] <. 1 



(A.6) 



There are multiple terms (or regime) to consider but the most important ones are when p 
is in the IR. For example, in the regime where p,q <^ k, the external momentum provides 
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the most stringent constraints and so we have that {t,ti,T2} ^ —1/k. So a typical term in 
Eq. ( [A .51 ) goes hke 

P Jo J -1/k ^1 J -1/k ^2 

where as promised the result of the time integral is to give 1/k^^^ with Va = 1 above. There 
are other terms of the same order as this one and others which are smaller (usually by powers 
of —kr). The integral over internal q will give k"^ such that the overall correct scaling of 1/k^ 
for (0^) is recovered while the p integral is infrared divergent. If we use a massive propagator 
instead of the massless one, we would get a - ~ ^ or in terms of A/^ a log divergence. So, 
barring cancelations the sunset diagram (Fig. P) scales like 

^-^HT^k,)^- (A.9) 



VZ^-V^3 



e 



where we include the external legs and the overall 6^ function. This two-loop diagram is 
therefore subdominant to the equivalent two-loop chain diagram by one power of 1/e. 
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